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SUMMARY 


The approximate formiala of a linearized solution for the pressure 
field generated hy a moderate rate of heat release is given. The anal- 
ogies between the pressure waves generated by heat release and those 
generated by (l) mass release, (2) piston motion, or (5) a two-dimensional 
body in a supersonic stream are established analytically. The exact 
solution of an idealized problem in which a finite amount of heat is 
released uniformly at a section of a tube with a given rate, large or 
small, is also constructed. Though this Idealized setup can be only 
approximately fulfilled in practice, the analysis does give an answer 
to a fundamental question; Given the rate of heat release at a section 
of a tube, how strong is the shock wave generated? A similar analysis 
is made for the pressure waves generated by a point source in three 
dimensions. Some applications of the theory are given. 


INTRODUCTION 


One of the basic problems in combustion aerodynamics is: What are 

the dynamic effects produced in a medixim as a result of heat release in 
the medium? In particular, how strong are the pressure waves generated 
by heat release and to what extent are they inqportant in a specific 
problem? Such questions arise naturally in the study of the spreading 
of autoignition, the transient development of detonation wave, ignition 
by compression, and many other time-dependent problems. The problem of 
estimating the pressure generated by heat release also occurs in other 
fields, for exeonple, in the study of spark discharge and thermally driven 
acoustical oscillation as well as of the mechanical effects ("blast”) 
produced in an atomic explosion (ref. l) . 

The basic mechanism by which pressure waves are produced by heat 
addition is slnq^ly this: When heat is added to a volume of gas, the 
density of the gas is, in general, reduced. This causes an expansion 
of the volume occupied by the heated gas. Ihe expansion of this volume 
produces the pressure waves . 
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To see this in a more quantitative manner, it is necessary to write 
down the fundamental equations governing the motion of the gas. Let 

p, p, T, and if he, respectively, the pressure, density, temperature 
velocity vector of the flow, all being seme functions of the position 

vector ^ and time t. let R and he the gas constant and specific 

heat at constant volume of the gas, hoth of which will he assumed constant 
in this preliminary study. The equations of continuity, momentum, and 
energy are 


1 + V . i? = 0 

p Dt 

(la) 

^ - VP = 0 

(lb) 

Dt p 




iflog ^ 

Dt'C pV T 

(Ic) 


vhere — represents the Stokes' derivative and q(^r,t^ is the rate 

of heat release per unit mass of the medium at r and at the instant t. 
When heat is added into the medium from an external source, q may be 
considered as given and equations (la), (lb), and (ic) together with the 
equation of state 

p = pRT (Id) 


form a system of four equations which govern the four unknowns p, p, 

m j 

T, and u. 


Wow it is observed that the rate of expansion of a given volume v 


of gas considered as a free body is measured by 


fj.t 


ds where s is 


the bounding surface of the volume and i? is the outward-drawn normal. 
Hence, by Gauss' theorem, it is also equal to 



u dv 


It then follows immediately from continuity equation (la) that whenever 

there is a change of density j there will be a change of the 

\P Dt y 

volume. Now, when heat is added to the medium, there will. In general. 
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"be a change of pressure and density in accordance with energy equa- 
tion (ic). (The relative amount of pressure change and density change 
resulting from heat addition can in principle always be determined in 
any specific problem by solving equations (l) jointly.) The change in 
density produces a change in volume occupied by the heated gas, which 
in turn generates the pressure waves. 


It. should be remarked that when heat is not added into the medium 
from external sources but is released by the fluid particles themselves 
the rate of heat release q must be considered as unknown and an addi- 
tional equation is required to describe its variation. In combustion, 
q is usually given as a function of the local temperature. In other 
cases, the rate of heat release of the fluid particles is specified. 

In any event, the system of equations (l) will then be considerably 
more complicated than it appears to be. 


Finally, one notes that the term Cy in equation (ic) is 

directly related to the entropy of the gas by the formula 


S - 



( 2 ) 


where the subscript o denotes some reference state and S is the 
entropy of the gas. 

For a moderate rate of heat release, the system of equations (l) 
can be linearized. The linearized theory is given in the following 
section where a "reduction theorem" is derived which enables one to 
reduce the problem of heat release in a tube to an equivalent problem 
of piston motion. An application of this reduction theorem to the 
one -dimensional case is given in the section entitled "APPIilCATION OF 
REDUCTION THEOREM. " When the rate of heat release is high, the 
linearized theory is no longer valid. In the section entitled "EXACT 
SOLUTION," an exact solution is obtained for the case when heat is 
released uniformly at a section of the tube at a constant rate (large 
or small) . The corresponding problem in three dimensions is given in 
the section entitled "THREE-DIMENSIONAL THEORY. " Seme applications of 
the theory are also given. 

The present investigation was conducted at the Department of Aero- 
nautics of The Johns Hopkins University under the sponsorship and with 
the financial assistance of the National Advisory Committee for Aero- 
nautics. The author is greatly indebted to Drs. Francis H. Clauser and 
Leslie S. G. Kovasznay for a few stimulating discussions and for their 
encouragement and to Dr. M. V. Morkovln for his interest and comments. 
The author would like to acknowledge his gratitude to Dr. Harold Mire Is 
of the NACA Lewis Fliglit Propulsion Laboratory who made several most 
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interesting suggestions. He is also grateful to Miss Ingeborg Busemann 
who carried out some of the calculations and to Miss Vivian O'Brien and 
Mrs. Sture Karlsson for their assistance in the preparation of the 
manuscript. 


SYMBOIS 


E 

H(x),H(t) 

L 

M 

mfr^t) 


cross-sectional area of a tube 
velocity of sound 

specific heats at constant volume and constant pressure, 
respectively 

energy per unit volume 

some arbitrary functions 

step function 
a characteristic length 
Mach number 

mass injected per unit volume i>er unit time 


n 


P 

Q 


Qo 


Q 

q(r,t) 
<l(x,y,z,t)[ 


R 


r 

r' 


normal to a control surface 
pressure 

rate of heat release, energy/sec 
constant rate of heat release (eqs. (46)) 
heating value of a mixture, energy/mass 

rate of heat release per unit mass, energy/mass/sec 

gas constant of a mlxtiore 
radial distance 

distance between field point and source point 
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r 

^1 

S 

St 

s 

T 

t 

U 

u 

-> 

u 

% 

^1 

Vs 

V 


x,y,z 


a 


radius vector 

radius of contact surface 

radius of spherical shock wave 

entropy 

transformation velocity (i.e., flame speed) 
hounding surface of volume v 
temperature 

temperature tending to “ in such manner that 
p#T*- = Constant 

time 

superficial velocity of a piston (see eq. ( 15 )) 
velocity (in one -dimensional case) 

velocity vector 

velocity of contact surface 

radial velocity (in three-dimensional case) 

velocity of flow immediately behind shock 

velocity of propagation of shock wave 

volume enclosed by control surface 

Cartesian coordinate of a point 

rate constant in parabolic heat-release law, equa- 
tion ( 47 )^ energy/sec3 


7 “ 
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6(x) 

5(x,y,z) 

e 

V 


p 

p* 


O) 

o)(y,z^t) 
Subscripts : 
c 
o 


1 


2 


increment of or change of 
&-f\inction in one dimension 
6-function in three dimensions 
any number , value of x 

vector normal to side of tube 
coordinates of source point 
density 

density tending to zero in such a manner that 
= Constant 

rate of heat release per imit area, energy/area/sec 
constant rate of heat release per unit area (eq. (20)) 
rate of heat release at x = |, energy/area/sec 

state at contact surface 

undisturbed state (i.e., that ahead of shockwave) 
except in case of variables representing rate of 
heat release to and Q (see list of main symbols) 

state immediately behind shock wave (which is also 
that immediately ahead of flame front in one- 
dimensional case) 

state immediately behind flame front 


LINEARIZED THEORY AND SOME ANAIOGIES 

Consider a medium in a uniform state with pressure Pq, density Pq, 
and temperature Tg. After the heat is added, the pressure, density, 
ten 5 >erature, and velocity induced can be written as 



MCA TN 5411 

P = Po 
P = Po 
T - 



T 

+ 6p 

(3a) 

+ 6p 

(5t) 

+ 5T 

(5c) 

bi? 

(3d) 


If the rate of heat release is not too hi^, — , — , and « 1 

Po Po Sq 

where Bq is the speed of sound in the undistxirhed medium. Equations (l) 
can then he linearized and 


+ V . 6u = 0 

St\Po/ 


— (hu) + = 0 

=,4-^ / o \7VoJ 


St 


S^-^p \ _ J_/5£\ = 
^^\rPoJ ^4\Poj CpTo 


5p ^ 5p _ ST 
Po Po ~ To 


Eliminate &p/Pq from equations (4a) and (4c): 

a |'5P.'\ + V . ei = ^ 


St\rPo/ 


^p^o 


(4a) 

(4h) 

(4c) 

(4d) 


(5) 


The equations governing the pressure and velocity fields are then given 
by equations (4b) and (5). Once the pressure field is known, the density 
and temperature field can be determined from equations (4c) and (4d) and 
the entropy spottiness from the linearized form of equation (2), 

^ ^ ^ (6) 

Cv Po Po 

It is now clear that "moderate rate of heat release" means — « 1 

CpTo^o 

where L is some characteristic length in the problem, for example, a 
length characterizing the dimension of the heating zone. 
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It is interesting to note that as far as the pressure aM velocity 
fields are concerned, the governing equations (eqs, (5) and (4b)) are of 
the same form - in the linearized approximation - as those which govern 
the pressure and velocity fields produced by injecting fluid (isen- 
tropically) into the medium. For, in the latter case, the momentum 
equation in the linearized form is the same as equation (4b) while the 
linearized continuity equation is given by 



(7) 


where m is the rate at which fluid (measured in mass/volume/second) 
is introduced into the medium. In fact, all one has to do to get the 
same pressure and velocity fields for the two cases is to match the 
parameter q/CpT^ with m/pg. This suggests that the noise generated 

in a nonuniform combustion (e.g., in turbulent combxistion) will radiate 
as a source field when examined at a large distance away from the com- 
bustion region, (it is not known if axay experiment has been carried 
out in the study of this phenomenon.) 


If &u is eliminated from equations (5) sud (4b) the pressure field 
is found to satisfy the wave equation 



It is clear from equation (8) that the pressure produced by heat addition 
depends upon the rate of heat release, and whenever there is a change of 
rate of heat release there will be pressure waves generated, a fact which 
one should have expected in the first place. 


Wow consider the gases contained in a tube. The x-axls will be 
chosen parallel to the axis of the tube. The heat added to the 
medlian q(x,y,z,t) is assumed to be given for t > 0. For t _ 0, 
q(x,y,z,t) is assumed to be identically zero. Now since 


r oo 

qU^y^z^t) 6(x - I) d| 


and equation (8) is linear (so that the superposition principle is valid), 
it is clearly sufficient to analyze the effect of heat release when 

q(x,y,z,t) = cs(y,z,t) 6(x - i), that is, when heat is added only at 

the section x = |. The quantity (s(y,z,t) has a physical meaning of 
its own. It is the heat generated per unit area per unit time at the 
section X = I, for 
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JJ caiy,z,t) 


dy dz = Total heat generated at 


X = I 


in unit time 


where the integration is taken over the cross section of the tube at 
X = |. Mathematically, this amounts to finding the elementary solution 
of equation (8), for once the pressure field produced by 0 )(y,z,t) at 
X = ^ is known, the general case of heat release by arbitrarily distrib- 
uted sources q(x,y,z,t) can be obtained by replacing the tenp o(y>z,t) 
in the formulas found for the pressxire and velocity field with 
S’^'i integrating them from | = -« to 

By a translation of the coordinate axes the plane at which heat is 
added can be made to be the plane x = 0. The pressure field produced 
by heat released at the plane x = 0 will then satisfy the differential 
equation 


82 

fbp \ 

_ 2 \ - 1. 

oj(y,z,t) 6(x) 

CVJ 

-P 

\7Voj 

" ° \7Poj 

Po^P^O 


and the boundary condition 

= 0 ( 11 ) 

avW 


at the vail of the tube, (Here v stand_s for the normal vector to the 
side of the tube.) Equation (ll), according to equation (4b), is an 
equivalent statement of the requirement that there should be no flow 
across the wall of the tube. If it is further assumed that initially 
at t = 0 there are no disturbances inside the tube, there are also 
the initial conditions 



(12a) 


(12b) 


The solution of the above problem can be facilitated with the aid of the 
"reduction theorem" given below. 
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Reduction Theorem 


Let heat be released at a section x = 0 of a tube of constant 
cross section and of infinite length at a rate of o)(y,z^t) units of 
energy per unit area i>er unit time. The pressure and velocity fields 
generated as a result of this heat release are identical to those pro- 
duced by two pistons at x = 0± moving away from each other (along the 

axis of the tube) with the same superficial velocity ^ ~ ^ ^ 

27 Po 


In this manner, the problem of pressure waves generated by heat 
release is reduced to a problem of piston motion. This reduction 
enables one to make use of the result of known theories (e.g., ref. 2) 
for the problem at hand. It is possible to generalize the theorem to 
a tube of, for example, finite length. This will be discussed later. 

The proof of the reduction theorem is very similar to that used in 
deriving the "impulse method" of solving a nonhomogene ous wave equation. 
(See, e.g., ref. 5>) 


Proof of Reduction Theorem 

First of all, observe that every solution of equation (lO) will be 
a solution of the homogeneous wave equation at all points except x = 0, 
since &(x) = 0 for x ^ 0. Next perform an integration of equation (lO) 
with respect to x from -e to e and then let e — >-0. If it is 

assumed that and are bounded near x = 0, 

^ t2 VtPo/ ay2 \T9o) \rPoJ 

a condition which can certainly be satisfied if o>(y,z,t) is a function 
smooth enou^ in y, z, and t, equation (lO) becomes 




o>(y,z,t) 


Po^^o 


(15) 


By symmetry, 6 p/tPq is an even function of x so that it is sufficient 
to consider the pressure field in the region x > 0 only. Also, by 

symmetry, (— = -(^ so that eqmtion ( 15 ) can be 

VSx n>oj^o_ V^x 7Po^=o^. 

written simply as 



1 d [7 - 1 g>(yjz>t) 
2,^2 stL 7 Po 


( 15 a) 
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In other words, if Sp/7Po is ^ solution of equation (lO), it must also 
satisfy condition (l3a) as x— >Of. Consequently, the solution 6p/7Po 
of equation (lO) satisfying conditions (ll) and (12) must satisfy the 
homogeneous wave equation 


/6P \ 
bt^V^o) 



= 0 


( 14 ) 


for X > 0 as well as conditions (ll)^ (12), and (l5a). But the 
pressure field in the region x > 0 produced hy a piston at x = 0+ 
moving with a superficial velocity 

(15) 

27 1^0 


in the x-direction will precisely satisfy equations (ll) , (12 ) , (l5a) 7 
and (l4) for x > 0.^ Moreover, since the pressure field generated hy 
a moving piston is uniquely determined in the linearized theory, it is 
also the only solution which satisfies equations (ll), (12) , (l5a) , and 
(l4). Consequently, the pressure field generated hy heat release and 
piston motion must he identical. The velocity produced hy heat release 
and hy piston motion are also identical for x > 0, since in hoth cases 
the velocity and pressure are related hy formula (4h) and 


StVTPo/ 


+ V 


5u = 0 


in this region. (Cf. eq. (5) and observe that q(x,y,z,t) = 

-i- u)(y,z,t) 6(x) and 6(x) = 0 for x ^ 0. Hence q = 0 for x ^ 0.) 
Po 

This completes the proof of the reduction theorem. 

It is clear from the above proof that the same conclusion will 
hold for tubes with nonuniform cross sections which are symmetrical 
about the plane x = 0. It also holds for tubes of uniform section of 
finite length provided that it is limited to the region ahead of the 
wave front reflected from the end of the tube. For the most general 
case, in which no such restrictions are in^josed, a similar but less 
simple analogy can be constructed. In this case the two pistons will 
be moving with different velocities, and they are so related that not 


^at it satisfies equation (l3a) follows from the momentum 
equation in the x-directlon (eq. (4b)), the boundary condition that 
there should be no flow crossing the face of the piston, and equa- 
tion (15) . 
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only is equation (l5) satisfied "but also the condition 



is fulfilled, (ihe equality of pressure is a 


direct consequence of momentum balance at the section x = 0 when all 
the second-order terms are neglected.) In the next section, applica- 
tion of this reduction theorem will he given. 


The reduction theorem was originally formulated as a means of 
reducing a problem of heat addition to a problem of piston motion. 
However, the fictitious pistons described in the theorem actually have 
a real physical significance. It has been shown in the Introduction 
that when heat conduction is neglected the flow field resulting from 
heat addition into a medium is really caused by the volumetric expansion 
of the hot gas resulting from the heating. In the case of heat addition 
at one plane (say, x = O) , the rate of this expansion can be calculated. 
In fact, the faces of the two fictitious pistons must be precisely the 
two interfaces which separate the cool gas from the heated gas , since 
they start from x = 0± at t = 0 and move out with a speed equal to 
the speed of the fluid particle. 

Dr. Harold Mirels pointed out that the reduction theorem states 
that the pressure and velocity fields produced by euidltion of tw units 
of energy per unit time per unit area are equivalent to those produced 
by two pistons at x = Ot moving away from each other with velocity 


6u 


7-1 CO 

2.J Po 


However, the rate of work done by such pistons is 2 Pq6u = 


7-1 


CD, 


whereas the heat input is co. Thus there are — cu units of energy 

7 

apparently unaccounted for. He went on to say that this is due to the 
step-function behavior of mass and energy at x = 0 (as can be deduced 
for the particular case discussed in the next section frcxti equations (22) 
to (24)). 


A more direct physical explanation of the energy that appears 
"unaccounted for" is the following. The cu units of heat energy per 
unit area per unit time added can be divided into two parts; That 
responsible for the setting up of the pressure and velocity fields 
(which, according to the reduction theorem, must be equal to 

2pq5u = ~ ^ co^, and that responsible for heating up the medium. On 

the other hand, in the case of the flow field prcaduced by the piston 
motion, all the energy goes into the setting up of the pressure and 
velociiy fields. The amount of energy that was previously unaccounted 
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for must then he the thermal energy stored in the hot gas- This can 
he seen from the following calculation. From the reduction theorem, 
it is concluded that the flow fields produced hy heating and piston 
motion are identical up to the faces of the pistons. In the case of 
heat addition, the hot gas filled the space between the faces of the 
two fictitious pistons. In a time interval 6t, the Increase in volume 
of the hot gas is 2ASu5t, where A is the cross-sectional area of the 
tube. The Internal energy stored in the hot gas is increased (during 
the same interval) hy 2A5u6t(pQ + 5p)C^(To + 5T) , that is, hy 


j ^ I A5u6t(Po + 6p) , 


or, to the order of the linearized approximation. 


simply hy j Po6u6tA. Hence the rate of increase in the internal 

energy stored in the hot gas per unit cross-sectional area at any instant 


t is 



— - — p /2_z_l iiL\ = which is exactly equal to 
7-1 2/ Po/ 7 


the amount of energy that appears unaccounted for if only the flow field 
produced hy the fictitious pistons is examined. (The above calculation 
is made for the case 03 = Constant across the cross-sectional area of 
the tube. For the more general case o) = o)(y,z,t), the same calculation 


applies except that one should write the integral 


dA instead 


of A in the calculation.) 


In the linearized theory presented here, only the effect of com- 
pressibility has been taken into account and the effect of heat conduc- 
tion has been neglected completely. Linearized theory taking heat 
conductivity into account was recently investigated by Wu (ref . 4) . 


APPLICATION OF EEDUCTION THEOREM 


As an application of the reduction theorem, consider the case 

where 


o)(y,z,t) = o>(t) 


(16) 


that is, the case in which heat is uniformly released at the plane x = 0. 
According to the reduction theorem, the pressure and velocity fields 
induced are the same as those generated by two pistons at x = Oi moving 
away from each other with velocities 


u(t) = ±l 


- 1 cp(t) 

27 p7 


( 17 ) 
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respectively. If the tube containing the gas is of constant cross 
section, the pressure and velocity fields produced by the pistons will 
be one dimensional, and they are given by the well-known formulas 


^ - 

7 - 1 


for 

X > 0 

(l8a) 

TPo 

27 

S'oPo 

sp ^ 

7 - 1 

^ k) 

for 

X < 0 

(l8b) 

TPo ~ 

27 

®-oPo 

Bu = 

7 - 1 


for 

X > 0 

( 19 a) 

27 

Po 

Su = 

y _ 

1 “(* ^ 

for 

X < 0 

( 19 b) 

27 

Po 


In particular, 
that is. 


if the rate of heat release is constant for t > 0, 


a)(t) = 



for t < 0 
for t > 0 


I 


( 20 ) 


or, written in terms of the step function, o)(t) = o:^(t) , the pressure 
field is given by 


5P _ 

7-1 «o H/'i 


for X > 0 

(21a) 

TPo 

27 aoPo \ 



5P ^ 

1 “o n(i 


for X < 0 

(21b) 

7Po 

27 SqPo \ 





This shows that, when heat is released at a constant rate of units 


of energy per unit area per unit time into the medium, two con^jression 
waves of equal strength are generated and propagate away from the heating 
zone. These compression waves have a strength, measiired in terms of the 
ratio of pressure jump 5p across the wave to the undisturbed pres- 


sure Pr,. equal to 7. ~ — It is thus seen that the strength of 

2 a^Po 

the pressure waves produced by heat addition to the medium is usually 
small. ^Thus, if ajt) = 10 Btu per square foot per second. 
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= 1,000 feet per second, and Pq = 2,000 pounds per square foot, 

then ^ = 0.2 X — ^ = 778 X 10-6 = O.O78 percent. ^ When the 

Po 1,000 X 2,000 ) 

rate of heat release is high, the pressxire waves generated may he quite 
intense. However, for a very hi^ rate of heat release the above formulas 
can no longer he used. In fact, the formula is correct only if the non- 

dimensional heat-release j)arameter « 1, since it is based on a 

®-oPo 

linearized theory. However, the siniplicity of the pressure and velocity 
field predicted from the linearized theory in the last instance leads 
one to suspect that an exact solution cam he found by replacing the 
infinitesimal pressure steps by two shock waves. Ihe main problem is 
again to determine the strength of the shock wave, given the rate of 
heat release o^j/a^Po. This problem will be considered in the next 
section. 


It is seen from equations (18) and (19) that the pressure is con- 
tinuous at X = 0, but the velocity is discontinuous there, since 




y - 1 g)(t) 
27 Po 


(S^^x=0- 


y - 1 g)(t) 
2 y Po 


The fact that the flow is leaving the section x = 0 at both x = CM- 
and x = 0- leads one to inquire if all the conservation theorems are 
really satisfied at x = 0. (it is certainly obvious that they are 
satisfied at all values of x ji^ 0.) 

It will now be shown that these laws are indeed satisfied at x = 0. 
For this purpose it is necessary to calculate the density and temper- 
ature field. Integrating equation (Uc) with respect to t. 
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By equation (4d) , 



Consider a control surface formed "by the planes x = e, x = -g (where 
e may he any number^ large or small ) f and the segment of wall of the 
tube between these planes. let A be the cross-sectional area of the 
tube. The conservation of mass and energy can be expressed as 



2p(e,t)u(e,t)A 


(24a) 



Substituting equations (3), (l8), and (19) into the above equations and 
neglecting the second-order small quantities (in accordance with the 
linearized approximation ) , it is found that equations (24) are indeed 
satisfied for all values of e. (By reason of symmetry of the flow field 

respect to x = 0^ the momentum equation formulated for the same 
control surface is automatically satisfied.) If € is taken as an 
arbitrarily small number, this establishes the validity of the conserva- 
tion laws at X = 0. If e is chosen large enough to enclose the whole 
disturbed field, it establishes the validity of the over-all conservation 
laws for the system. Note that these conservation laws will not be 
satisfied if the pressure and velocity fields are related to the rate of 
heat release a> in any other manner than by equations (l8) and (19) • 

There are two aspects of this analysis which deserve some criticism. 
First of all, it is clear from equation (25) that, in general, 6 T/Tq 
is infinite at x = 0. This is, in fact, a direct consequence of the 
addition of a finite amount of heat in a plane. Since a "plane of gas" 
(instead of a volume) possesses no heat capacity, the temperature must 
become infinite. Now, if 5 T/Tq is infinite, the linearization must 
break down at x = 0 so that the solution really does not represent any 
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physical state of affairs at x = 0. Secondly, the density 6p/Po^ 
according to equation (22) , becomes at x = 0. Here there is not 

only an Infinity hut also a negative density which is not even physi- 
cally conceivable. It will be seen in the next section that the nega- 
tive density is really a consequence of linearization. 

2 

Despite the absurdity of the behavior of the solution at x = 0, 
all the other conclusions derived above, such as the dependence of the 
strength of the pressure vave on the rate of heat release (eq. (2l)) are 
to be trusted because only this depen(?.ence will insure the conservation 
of energy, mass, and momentum locallv at all points as veil as for the 
vhole system (i.e*, over-all balance}* This state of affairs vill be 
clarified when an exact solution is constructed in the next section- 


Other evidence which indirectly justifies the value of this solution 
is found when use is made of the above solution to get the pressure and 
velocity fields produced by heat sources q(x,t) which are not distrib- 
uted in one plane but over a region (volume) of the tube- According to the 
rule given (see eq« (9) and the paragraph following it) , it is necessary 
only to replace x in formulas (l8) , ( 19 ); (22), and (25) by x - | 

(i.e., perform a translation) and then replace the function o) in these 
formulas by Poq(l^t) d|, and, finally, integrate i from -00 to 

Thus, the pressure field is given by 


^ 7 - ^ 
TPo 27 




But, since q(^>t) = 0 for t ^ 0, the last formula can also be written 
as 


8p 

TPo 




CpTo 


d| + 


^x+a^t q(|,t + 


d| 


(25) 


which is well 
equation 


known as the correct solution of the one— dimensional-wave 


5t2\7Po/ 


a 2 /SP \ = A 

q(x,t) 

° dx^VTPo/ 

1 

0 

1 


subjected to the initial conditions 



^To be STore, to begin with, it is physically impossible to add a 
finite amount of heat in a plane (instead of a volume) of gas. However 
this is beside the point because this state of affairs can at least be 
approximated in practice- 
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Likewise, the velocity, density, and temperature fields are given hy 

, rx+aot qf|,t + 

(26) 

pXo ^ X ^’p-*^o 


&U = i 


X a(S.t - ^ 1 r*«ot ^ 

2 „ 4 - Cr,To ” 2 


x-a^t 


5 p _ 1 


Po 2aQ 


r^X 


x-a^t 


i(g.t - 


CpTo 

x+a^t q(|,t + 


d| + 


d| 


q(x£^ (27) 

0 *^"^0 


OT 

To 


- 7-1 


2a.r 


^x q 


x-aQt 


- — 7n ^ — ^ + 

CpTQ 


>x+aot q[|,t + 


CpTo 


d| 


^ dt 

0 


(28) 


Note that in the expressions for the temperatiore (eq. ( 28 )) and density 
(eq. (27)) the "infinities," which were originally contained in expres- 
sions (22) and (25), disappear after superposition. 


To get some idea as to the pressure and velocity distributions inside 
a narrow (but nonzero) band of heating zone, a calculation is made for 
the case 


q(x,t) 


0 

lo 


for t < 0; also for t > 0 if jx| > g 
for t > 0 and jx| < g 


where e is some small number. Typical pressure and velocity distribu- 
tions for various times are shown in figure 1. 


The solid lines in figure 1 are the pressiare and velocity distribu- 
tions. The dotted lines are the construction lines. Points between aa' 
do not realize the heating zone is cut off at x = ±e and think that it 
extends from -» to ». By symmetry, therefore, the fluid particle tends 
to move toward neither the left nor the right. Hence in this region 
6u = 0, and the density, by the continuity equation, remains unchanged, 
while the pressure increases linearly with the temperature (constant- 
volume combustion) . Points between ab do realize that the heating zone 
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is cut off at X = € but do not know if it is also cut off at x = -€. 
The pressure and velocity differ at different points in ab depending 
on how much each point knows about the conditions at each side of the 
point X = €• The same remarks hold for points between a’b*. 


Before concluding the discussion of the reduction theorem, one 
should mention that when heat is released Txniformly at one plane, that 
is, o)(y,z,t) = 03 (t), it is possible to construct a third analogy, 
namely, that between the pressure waves generated by heat release and a 
two-dimensional airfoil in supersonic flow. 5 For the reduction theorem 
shows that the pressure field produced by heat release satisfies the 
system of equations (l4), (15), (l2a), and (l2b), which in the present 
case become 


- a 2 = 0 


^ /bp \ 

- 1 5 

7 - 


x=0t 

o * 

CVJ 

1 


(29a) 

(29b) 


f— 'I '0 

( 29 c) 

- 0 

atlrPo; ^0 

( 29 d) 


On the other hand, if a symmetrical two-dimensional airfoil y = ±f( 
in a uniform supersonic stream of Mach number M is considered, the 
pressure field produced by the airfoil satisfies the following set of 
equations : 


a2 /5p \ _ 1 ^2 \ ^ Q 

^x2\^o/ Mq 2 _ 1 dy2\XPo/ 


(50a) 


a /5P 

SyV^o 


y=ot 


= T A f(x) 


( 30 b) 


^The author is indebted to Dr. L. S. G. Kovasznay for suggesting 
this analogy in the early stage of this analysis. 
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( 30 c) 


h /^V \ 

SxVtPo/ 


x=0 


0 


( 30 d) 


(According to the momentum equation in the y-direction eq. (50h) is 
equivalent to the statement that the flow at the surface of the airfoil 
must he tangent to the airfoil. The last two conditions are equivalent 
to the statement that there should he no incoming waves in the flow 
field.) 

A comparison of equations (29) and (50) shows that there exists an 
analogy between the pressure field developed hy the uniform heat addi- 
tion at a plane and the two-dimensional symmetrical airfoil in a super- 
sonic stream. In fact, it is necessary only to imagine the time coordi- 
nate as the space coordinate in the direction of the supersonic stream 

and choose IV^ = 1 + and the shape of the airfoil according to the 

ao^ 

formula f(x) = ^ ~ ^ i ^ . 

27 1 + a.2 Po 


EXACT SOLUTION 


Let heat he uniformly released at the section x = 0 of a tube of 
constant cross section and of infinite length at a constant rate of 
units of energy per unit area per unit time, (Hence^ if the total heat 
released per second is Q and the cross-sectional area of the tube is 
= q/a.) When the nondimensional heat-release parameter co^/a^pQ vas 

y — 1 ^o 

small compared with unity, two con^jresslon waves of strength ^ — a^^^ 

propagating away from each other with the velocity of sound a^ of the 

xmdlsturhed medium, were found to he generated at x = 0 (see eq. (21)) . 
When <%/8-oPo small, equation (2l) is no longer valid. Further- 


more, it is expected that, instead of the two infinitesimal pressure 
steps, two shock waves would he generated. The question is again to 
determine the strength of the shock wave generated, given the rate of 
heat release per unit area 
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The linearized solution (21 ) , though invalid for a hi^ rate of 
heat release, does siiggest the following propositions: 

(1) The shock waves generated at x = 0 are of equal strength and 
will maintain their strength as they propagate away from each other. 

(2) The pressure between the two shock waves is uniform and equal 
to the pressure immediately "behind the shocks. 

(5) The velocity field "behind the shock waves is discontinuous at 
X = 0 and is an odd function of x. It is uniform for the regions x > 0 
and x < 0, respectively, and assianes the value of the velocity of the 
gas immediately "behind the shock in the region x > 0 and x < 0, 
respectively. 

In addition to these propositions, due account will "be taken of the 
motion of the fluid particles which is neglected in the linearized theory. 
It would "be expected that: 

( 4 ) Two contact surfaces which form the "boimdaries separating the 
hot and cool gas are generated at x = 0 at the instant when heat is 
first added to the medium and move away from each other with the velocity 
of the fluid particles. 

(5) Since the hot gas originated from a single plane x = 0, the 
temperature of the gas between the two contact surfaces must tend to 
infinity. 

(6) For the same reason, the density of the hot gas between the two 
contact surfaces must tend to zero in such a manner that pT = Constant 
in accordance with the gas law pT = p/R and proposition (2) . 

Wiese propositions can be substantiated by a formal argument based 
on dimensional reasoning. The argument goes as follows: The undisturbed 

medium can be characterized by two of its thermodynamic state parameters, 
say, the pressure Pq and temperature Tq. Since the velocity of 

sound a^ in the undisturbed medium is uniquely related to the tenqier- 
ature Tq, Pq and Uq will be used as the two parameters characterizing 

the undisturbed medium. The strength of the shock wave can be described 
in terms of the pressure ratio Pi/Pq across the shock, where p^ is 

the pressure immediately behind the shock. It is clear that, in general, 
the strength of the shock wave depends on the rate of heat release per 
unit area the state of the undisturbed medium being characterized 

by Pq and aQ as well as by the time t. That is, 

= F(c%,aQ,PQ,t) 


( 51 a) 



22 


MCA TN 5^11 


Since the viscous and heat-conductive effects have been neglected in 
this preliminary study, these variables do not enter into equation (31a) • 
Also, for uniform heat release in a tube of constant cross section the 
flow field is one dimensional so that the dimension of the tube does not 
enter into the problem as a relevant characteristic length. Wow equa- 
tion (51a) must be dimensionally correct. But the four variables o:^, 

Pq, Uq, and t can only be combined into a single nondimens lonal param- 
eter, namely, ^ , which does not contain the variable t. Consequently, 
®-oPo 



(31b) 


that is, the shock strength must be Independent of the time t, which 
proves proposition (l) . Clearly this conclusion is actually a direct 
consequence of the fact that there is neither a characteristic time nor 
a relevant characteristic length in the problem. 


Likewise, the pressure, density, ten 5 >erature, and velocity of the 
flow behind the shock waves will be fimctions of cq^, a.^, Pq, and t. 

In addition, they can be fimctions of the position x. The five vari- 
ables cut,, Uq, Pq, X, and t can be ccatibined to give two independent 

nondimensional parameters, namely, and _2^. Therefore, 

aQt 


and so forth. 


p(x,t) 


( 32 a) 

Po 

a-oV 

u(x.t) 

1 

N 

PJH 

II 

( 52 b) 

^■0 

^V^oPo 


words • 

the flow field must be "conical. " 

Intro- 


ducing a new independent variable ^ it is a sin 5 )le matter to 

reduce the governing partial differential equations (l.e., the continuity, 
momentum, and energy equations) to a system of ordinary differential 
equations with the independent variable It turns out that this 

system can be integrated without difficulty so that explicit solutions 
are obtained for the pressure, velocity, temperature, and density fields. 
The constants of integration are to be determined from the boundary 
conditions at the shock. However, a simpler approach is to assume a 
possible flow field which embodies all the features in propositions (1) 
to (6) and then to verify that it indeed satisfies all the conservation 
laws at all points . 
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For the purpose of formulation, suppose for the time being that the 
strength of the shock waves is known. Then the pressure P]^, velocity u^^, 
tenqjerature Tp, and density immediately behind the shocks are also 

known. Since the shock strength is invariant with time, pp, up, pp, 
and Tp will also be independent of time and are therefore constants. 

To satisfy proposition (2), assume 


p(x,t) 

= Pi 

for 

|x| < Vgt 

where Vg is the shock speed. 

Proposition ( 5 ) states ' 

u(x^t) =- 

fup 

for 

0 < X < Vgt 

i 

-Ul 

for 

-Vgt < X < 0 


(55a) 


(55t) 


As a cOTisequence of this formula and proposition (4) , the two contact 
surfaces which separate the hot and cool gas must be traveling with the 
velocity tup, respectivejiy. (See fig. 2.) Proposition ( 5 ) then 
requires that 


where T*<^ 
sition (6) 


T(x,t) 


Ti 

for 

T* 

for 


uit < |x| < Vgt 
1x1 < upt 


V II J 

is to be taken as a tenperature tending to infinity, 
requires that 


p(x,t) 



for Upt < jx| < Vgt 
for |x| < Upt 


(55c) 


Propo- 


(55d) 


where is to be taken as a density tending to zero in such a manner 

that pi«T* = 

R 


It is clear that this solution satisfies all the conservation laws 
for any shock strength at any point in iiie flow field except perhaps at 
X = 0. For, if the flow field in the region x > 0 at any instant t 
is considered, it consists of a shock wave at x = Vgt, followed by a 
flow with uniform pressure and velocity consistent with the shock rela- 
tions, and a contact surface at x = upt, which moves with the fluid 
velocity. Such a flow field certainly satisfies all the conservation 
laws in the region x > 0. Hence, what must be examined is the question 
of whether the conservation theorems are also satisfied at x = 0. If 
they are satisfied there, then the conservation theorems for the system 
as a whole (i.e., over-all balance) will also be satisfied and the 
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converse is also true. It is precisely this consideration that deter- 
mines the strength of the shock, vave. 

Take a control surface consisting of the planes x = e, x = -e, 
and the segment of the wall of the tube between them. The conclusion 
will he the same whether e is chosen to he an arhitrarily small number 
or a number large enough to Include the entire disturbed field. It is 
found that the equations representing the conservation of mass and 
momentum are always identically satisfied no matter what is the shock 
strength, while the energy equation is satisfied if and only if the 
strength of the shock waves is such that the following condition 

pm = 

is fulfilled. These statements can be proved easily if e is chosen to 
be less than Vgt. When c is taken to be larger than Vgt the calcu- 
lation is more Involved, since use has to be made of the shock relations. 

It will be shown how the calculation proceeds for this more Involved 
case (which incidently establishes the conservation laws for the system 
as a whole) . 

Take e to be some fixed number greater than Vgt. Then at time t 
the shock waves will be at x = Vgt and the contact surfaces, at x = u^t 
(see fig. 5) . The total mass inside the control surface consists of three 
parts: That between the two contact surfaces 2uitAp*, which is negligible 

since p* tends to zero eventually, that between the contact surface and 
the shock wave preceding it 2(Vgt - uqt^Apq, and that between the shock 
wave and the boundary of the control sxirface 2(e - Vgt)Apo. Thus, the 
total mass inside the control surface is 

2uj^Atp* + 2(Vgt - uit)Api + 2(e - Vgt)ApQ = 2u]^Atp* + 2eApo 
where use has been made of the shock relation 


Pl(Vs - ui) = PoVg (55) 

(i.e., the continuity equation at the shock). Now, since p* is actually 
zero, the total mass contained in the control surface is Invariant with 
time for all values of t and is eqixal to 2eAp^, the mass contained 

Inside the control surface before heat is released to the medium. This 
demonstrates that the mass is conserved. 
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The momen'tum equation, being a vector equation, is automatically 
satisfied by virtue of the symmetry of flow with respect to the 
plane x = 0. Then consider the energy balance. 

At the instant t, the energy inside the control surface also 
consists of three parts. It is given by 

2uitAp*(CvI^ + I ui^) + 2(Vgt - Uit)Api(CvTi + | u^^) + 


2 (e - Vst)ApoC^To = PlU^^tA + 2 €ApoC^Tq + 


2tAVgPo(c^Ti + I ui2 _ c^To) 


since p*- tends to zero in such a manner that p*T* = Pi/R. Again use 
has been made of the continuity equation at the shock (55) . The 
term VgPo(c^Ti + i uq2 - C^To^ can be converted into p^u^ if one 
makes use of the shock relations 


Cp(Tl - To) = VgUi - I ui2 


( 56 ) 


Pi - Po = Po''^s^l 


( 57 ) 


1 



Vs = ^1 


which are nothing but the energy, mcanentum, and continuity equations at 
the shock. Thus, the total energy contained inside the control surface 
at the instant t is given by 


— ^ PpUptA + 2eApoC^To + 2tApiUp = 2eApoC^To + ppUptA 

y - 1 ' 

The increase of energy in the control surface from the instant t to 
the instant t + &t is therefore (27/7 - l)piUpA5t where 6t is 

chosen so small that Vg(t + 5t) is still less than e. But the total 
amount of heat released in time t is tUpAbt. Consequently, the energy 
will he conserved if and only if 
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Pl^l = ^ 


which proves equation (5^) . 


One is now in a position to relate the shock strength w”ith the rate 
of heat released per unit area cju^. It is well known that the drift 
velocity U]^ produced by a shock wave of strength Pi/Pq 


^1 



a 


o 


zJlA El + 7- 1 

2/ Po 27 


( 58 ) 


Substituting equation ( 58 ) into equation (5^)^ the following relation is 
found between the strength of the shock wave generated and the parameter 
of rate of heat release o^^/a-oPo- 


03, 




2 £i/h 

7-1 pMPo 




^oPo 


/ y + 1 Pl~ 
V 27 Po 


- 1 


27 


(59) 


One can then solve the quartic equation for Pq/Po terms of (%/aQPQ. 
A plot of Pi/Po versus ®q/^oPo given in figure 4. Obviously, this 
curve has been constructed by finding the values of (JC^/aoPo ^ 

sequence of values of Pi/Po- results of this calculation are pre- 

sented in the following table: 


Pi 

Po 

a>o 

®'oPo 

Pi 

Po 

% 

®oPo 

2 

7.54 

10 

152.5 

4 

51.8 

15 

290.9 

6 

65.2 

50 

1,867 

8 

105.7 




It is interesting that for a very hi^ rate of heat release the strength 
of the shock wave varies with cq^/^'oPo according to a two- thirds law: 


Po 


0.525 



2/5 


(40) 
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(The coefficient 0.525 corresponds to a value of y of 1.4*) Once 
Pl/Po is known as a function of %/s^o^^o^ complete flow field is 
defined and given by equation (55) • 

Application of the foregoing results to the approximate estimation 
of pressiire waves generated by autoignition of a mixture in a tube is 
discussed in the section entitled "SOME APPLICATIONS OF THEORY. " 

THREE-DIMENSIONAL THEORY 


In the present section^ the linearized theory of pressure waves 
generated by heat addition in three dimensions will first be discussed 
briefly. Then the construction of an exact solution in three dimensions 
which may be useful in predicting the asymptotic strength of the shock 
wave generated by a closed flame front expanding uniformly will be 
discussed. 


According to the linearized theory^ the differential equation 
governing the pressure field due to a moderate rate of heat release is 
given, by equation (8) , The solution of equation (8) in an open space 
is well known: 


Sp(x,y,z,t) 

7Po 




f 1 a 


1 

8 

1 

CU| 
ct 1 

S'^o 



d| dT] d^ 


( 1 ^ 1 ) 


where r' = \/(x - + (y - il)2 + (z - is the distance between the 

field point (x,y,z) and the soiirce point 


The special case in which heat is released at a given point in 
space, say the origin of the coordinate system, is of particular interest. 
This means that 


q(x,y,z,t) = Q(t)6(x,y,z) 


( 42 ) 


where 6(x^y,z) is the 5-function and can he considered as the limit, 
as e — >0, of the function 


0 


for 


x^ + y^ + z^ > 


&(x,y,z) 



for x^ + y^ + z^ < 


( 45 ) 
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The function Q(t) has the physical significance of heing the heat 
release per unit time. Substituting equations (1+5) into equation (l+l) 
and evaluating the integral, 

5p(x,y,z,t) _ 1 1 ^ 

r 8t 


where r is the distance vx^ + 

A comparison of equations (21) and (1+1+) shows that one of the 
essential differences between the one -dimensional and three-dimensional 
cases is that the pressure waves generated in the former case depend on 
the rate of heat release (per unit area) while those generated in the 
latter case depend upon the time rate of change of rate of heat release. 

In practice, it is of some interest to know the pressure and 
velocity fields produced by the sudden addition of a finite amount of 
heat into the medium. Thus, the flow field produced by a spark discharge 
is of this nature, althou^ actually a linearized theory will not be 
adequate to describe this phenomenon with precision. Assume the varia- 
tion of rate of heat release with time as shown in figure 5- I++ other 

words, at t = Of, the rate of heat release Q(t) increases suddenly 
from zero to a very high value and then decreases again to zero in a 
short interval of time. The pressure waves generated, according to 
equation (41+) , will vaiy with the derivative of Q(t) and will there- 
fore consist of a very steep compression front followed immediately by 
an expansion-compression zone (fig. 6) . The velocity field produced 
may be calculated from equation (5) • The distribution of velocity field 
along any radial line is 

u. = i i + _k. i A 

4rt CpTo 4itao r ^t 





CpTo 


(44) 


Thus, the velocity distribution in the immediate nei^borhood of the 
origin behaves like an incompressible source field. Moreover, the 
radial distribution of the velocity there varies like Q (near field in 
fig. 6), while that at a large distance away bears the same relationship 
with the pressure as that existing in the theory of plane wave (far field 
in fig. 6) . 

Next examine the possibility of constructing an exact solution in 
three dimensions as was done before in one dimension. First, consider 
the case in which heat is released at a uniform rate of Qq units of 
energy per second at the origin. As in the preceding section, the 
undisturbed medium can be characterized by its pressure p^ and velocity 
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of sound Hq. Then the strength of the shock wave generated, measured 
in terms of the pressure ratio Pi/Po across the shock wave, will, in 
general, he a function of Qo^ ^-o-* Po-> time t. That is. 


= F(Qo,ao,Po,t) (46a) 

Now the four variables Qq, Qq, Pq, and t can only he combined into 
a single nondimensional parameter Since equation (46a) must 

he dimensionally correct, one has 


Po 



(46h) 


Consequently, in this case the shock strength must he a function of 
time - a fact which complicates greatly the construction of an exact 
solution. Since it is physically apparent that increasing the rate of 
heat release will increase the strength of the shock wave generated, it 
is concluded from equation (46h) that the shock wave generated must 
decay with time. A little reflection reveals immediately that the basic 
reason why the shock strength should depend on t is the existence of a 

characteristic time . / — in the problem. It is also clear imnediately 

V ao^Po 

that in order to produce a shock wave whose strength is invariant with 
time one must add heat to the medium according to the law 

Q(t) = at2 (47) 

where a now has the dimension energy/time^ . For, in this case, the 
rate at which heat is released will he characterized by a and the shock 
strength Pi/Po should he a function of a, Bq, Po> and t instead. 

But the four variables a, aQ, Pq, and t can only he ccmhlned into a 
single nondimensional parameter a^o^Po which does not contain t. 
Consequently, the shock strength Pi/Po will he a function of 
hut not of t; that is. 



The construction of an exact solution will he attempted for this case. 

The "parabolic law" (equation (4?)) is actually of some practical interest. 
Thus, if a spherical flame propagates into the fresh gas with a constant 
speed, the rate of heat generated by combustion is proportional to t . 
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Moreover, if the flame speed is small, the heat may he thou^t of as 
being released at the center of the sphere. 

For the present case, the pressure, velocity, density, and temper- 
ature field must be a function of a, aQ, Pq, r, and t. Just as in 

the derivation of equation ( 52 ), a dimensional reasoning leads to the 
result 


Po 



y 



(49a) 



(49b) 


and so forth, so that the flow field is again conical. It is also known 
that there must be a spherical contact surface behind the shock wave 
which separates the heated gas (i.e., the "fireball") from the cool gas 
surrounding it. If the position of this spherical contact siirface at 
any instant t is denoted by its radius then, in general, 

r^, = C(a,ao,pQ,t^ . In nondimens ioml form, this equation can be 
written as 



so that the contact surface must move out with a uniform velocity, or 
the fireball must expand at a uniform rate. Consequently, the flow field 
outside the contact surface will be exactly the same as that produced by 

a solid sphere expanding at a uniform rate of a^Gf — v — \ But the air 

V^o^Po/ 

waves generated by a uniformly expanding sphere have been solved by 
Taylor (ref. 5)^ so that the shock strength will be determined as soon 
as the rate of expansion of the fireball is related with the nondi- 
mensional heat-release parameter a/ao^Po- find this relation it is 

necessary to construct first a solution which is valid inside the 
fireball. 

Except at the origin r = 0 where heat is released, the fundamental 
hydrodynamic equations governing the flow Inside the fireball are 


^In fact, the dependence of pressure and velocity field on r and 
t in Taylor's solution satisfies precisely requirement (49). 
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St 


+ 


L 1 




0 


( 51 a) 



(511) 


^ P(^VT . i 


+ 


i_ i_ 

1*2 Sr 


r^pup 




0 


( 51 c) 


p = pRT 


( 51 d) 


Just as In the preceding section, it is expected that the temperature T 
in the fireball tends to infinity, while the density p of the gas there 
tends to zero in such a manner that 


pT = p/R (52) 

If this assumption is made, then continuity equation (51a) is automati- 
cally satisfied provided that Ur is finite. The momentum equation (51b) 
will then he satisfied if p is Independent of r. But according to 
equation (49a) , if p is independent of r it must also be independent 
of t. Hence, it is assumed that 


P = Pc 

where p^ is the pressure at the contact surface. Finally, since 

p P 

pT = the energy equation will be satisfied provided that r^Uj. is 
R 

a function of t only. To satisfy equation (49b), it is necessary that 
Ur = (Constant) t^/r2. in fact, if Uc is the velocity of the contact 
surface. 



because, at r = r^ (= U(.t) , Uj. must be equal to u^.. (Note that Ur 
is Indeed finite for all values of r > 0.) By symmetry, Uj. = 0 at 
r = 0. 

Summarizing the results, the flow field inside the contact surface 
is given by 


P = Pc 


(53a) 
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u 


c 



for r > 0 




=r 0 at r = 0 


T = 

P = P* 


( 5513 ) 

(55c) 

(55d) 


where T* tends to <» 


and p* tends to zero in such a manner than 


pm = Pc/R 


(5^) 


Note that inside the fireball the velocity distribution is like an inccan- 
pressible source field (just as that indicated by the linearized solution, 
eq. (45))- 

Up to now, the conservation laws have been satisfied at all points 
inside and outside of the fireball, except at the point r = 0 itself. 

An examination of the conservation laws at this point will enable one to 
relate the velocity of expansion of the fireball u^, and the pressure 
in the fireball Pc with the heat-release parameter a (see eq. (55))- 
If Uq is assumed to be known for the time being, Taylor's solution will 
then give the shock strength P 3 _^Pq and the flow field outside the fire- 
ball as well as the value of p^,. It will be shown that the correctly 
assumed u^. must be that which yields a Pc consistent with equa- 
tion (55)» 

Consider a spherical control surface of radius e about the origin. 
Take e so small that the control surface lies completely inside the 
fireball. The continuity and energy equations formulated for this control 
surface are: 


^t 



p*4rtr^ dr + 



= 0 



i Up^Vrtr^ dr + 


/ 




cm + 


+ Pc"^c(^) 



NACA TN 5411 


35 


The first equation is automatically satisfied since p* ultimately tends 
to zero. The second equation can he rewritten as 

so that it will he satisfied if aud only if 

PcUc5 = a (55) 

7-1 

Finally^ the momentum equation is automatically satisfied hy virtue of 
symmetry. Consequently, all the conservation laws are satisfied every- 
where provided that the flow field outside the fireball is described hy 
Taylor's solution, that inside the fireball, hy equation (55) ^ and equa- 
tion (55) is satisfied. It is often more convenient to use the nondi- 

mensional form — 2 — = of equation (55)- 

ao^Po 7 - 

The table below is essentially a reproduction (in the notation of 
the present paper) of table I of Taylor's paper (ref. 5) ^ except that one 
more column ( the last column) has been added to give the corresponding 
values of ay^Lo^Po in accordance with equation (55)* la the table Mq 
is the I^ch number of the contact surface and r]^ is the radius of the 
shock wave. Other notations have already been Introduced earlier. 


© 

© 

® 



© 

© 

Me 

^1 

^c 

Pi 

Po 

^c 

^0 

Pi 

Pc 

a 

Pc 

Po 

ao^Po 

0 

.2 

.4 

• 5 
.6 

• 7 
.8 

1.0 

1.2 

1.4 

1.6 

1.8 

2.1 






0 

.592 

5.89 

8.75 

17.74 

54.7 

64.9 

212 

651 

2,050 

7,250 

40,000 

00 

4.95 

2.44 

1.950 

1.765 

1.505 

1.392 

1.256 

1.182 

1.155 

1.105 

1.085 

1.060 

1.000 

1.005 
1.050 
1.169 
1.565 
1.629 

2.400 i 

5.59 

5.60 

9.06 

17.95 

00 

0.205 

.410 

.525 

.658 

.761 

.891 

1.180 

1.520 

1.955 

2.560 

5.598 

00 

^ 

0.928 

-775 

.750 

.749 

.755 

.774 

.811 

.847 

.887 

.917 

.92 

.93 

1.075 

1.295 

1.400 

1.569 

1.808 

2.105 

2.959 

4.250 

6.52 

9.89 

19.7 

00 
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Column (j^ of the table gives the heat-release parameter, while 
column (5) gives the strength of the shock wave produced. The dependence 
of the shock strength on the heat-release parameter is also 

shown in figure 7* Column (4) of the table gives the velocity of expan- 
sion of the fireball relative to the sound speed in the undisturbed 
medium. Column gives the pressure inside the fireball. Column 
gives the ratio of the radii of the spherical shock and the fireball. 


SOME APPLICATIONS OF THEORY 


Consider first some possible applications of equation (59) • It is 
clear that in practice there are cases in which heat is released at a 
uniform rate in a limited region (or a narrow band) of a tube. If the 
axial dimension of this region is small compared with the length of the 
tube, equation (59) would be expected to give approximately the correct 
strength of the shock wave produced by the heat released. If the tube 
is of infinite length then equation (59) would be expected to give the 
asymptotic value of the strength of the shock wave developed, whatever 
is the size of the region at which heat is actually released (provided 
that this region is finite) . A similar idea applies to some ignition 
problems. Consider a tube of infinite length containing some combustible 
mixture. Suppose that at t = 0, for one cause or another, ignition 
begins at one section, say at x = 0 , of the tube. Thus, the mixture 
may have been ignited by a hot surface or a grid, or automatically because 
of the existence of a local high- temperature region. A flame is developed 
and the combustion will tend to spread out into the fresh gas (see fig. 8). 
If the flame spreads out at a uniform rate, the total amount of heat 
generated as the result of combustion will be linearly proportional to 
time t and the rate of heat release is therefore constant. Furthermore, 
if the flame spreads out at not too high a speed, compared with the 
velocity of sound in the surrounding medium, the strength of the shock 
wave generated would be expected to be approximately the same as if heat 
had been released at a constant rate at the section x = 0 . In fact, one 
can calculate the strength of the shock waves generated. Thus, supposing 
that the two flame fronts propagate away from the ignition plane x = 0 
with a constant speed S-^., (the transformation velocity), the total amount 

of heat generated in time t is 2pQ_S-^AQt where is the density of 

the medium behind the shock waves, A is the cross-sectional area of the 
tube, and Q is the heating value of the mixture (in energy per unit 
mass) . It follows therefore that the rate of heat release per unit area 
is constant and given by 


0^ = 2piS^Q 


(56) 
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and the strength of the shock wave generated can he solved from 
equation (59) 


_2_ !if£i . i\ 

7-1 PqVPo / ^ 2 Pl^tQ 
V + 1 Pi ^ 7 - i ®’oPo 
27 Po 7+1 


( 57 ) 


where is related to Pi/Pq Rankine-Hugonlot relation 


£l = 

Po 


1 + 


7+1 £l 

7 - 1 Pq 


7+1 + ^1 

7 - 1 P^ 


Naturally, one likes to know if this gives the correct ansver. Fortu- 
nately, in the present case there is an independent method of calculating 
the strength of the shock waves generated. Note that since the flow field 
must he symmetric about the plane x — 0, the same flow field would have 
been produced if the tube had been closed off at x = 0 (fig* 9) • Now 
a flame propagating away from the end of a tube must generate a shock 
wave of such a strength that the boimdary condition at x = 0 is satis- 
fied. Let the shock wave have the strength Pi/Po >^here p^_ is the 
pressure immediately behind the shock. The shock will then induce a 
flow giving rise to a drift velocity (fig. lO) 





7^1 Pi 7 - 
27 Po ^ 27 


1 


(58) 


where blq is the velocity of sound in the medium ahead of the shock 
wave. Since the flame is assumed to be propagating with a constant 
speed S-(^ relative to the mediimi, it will be seen to propagate with 

an apparent speed of U]_ -h S-^. Now the flame itself induces a flow 
behind it. From the viewpoint of an observer riding on the flame, the 

burned gas is leaving the flame with a velocity equal to — where 

is the density of the gas ahead of the flame (i.e., that behind the 
shock wave) and p2 is that behind the flame. Hence, from the viewpoint 

of an observer in the laboratory (i.e., one who is fixed with respect to 
the undisturbed medium ahead of the shock waves) , the burned gas will be 
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moving away from the wall 


X = 0 with a speed equal to 


ui + St - St 


But the boundary condition at the wall dictates that the velocity of the 
flow at the fixed wall is zero. This condition is satisfied if 
Pi 

ui + St - St -^ = 0 or 


“1 = - l) (59) 

which together with eqixation (59) determines the shock strength since 
Q±/P 2 ^ fixed ratio once the mixture is specified, Now^ if u^_ is 

small compared with the velocity of sound aj in the mediiom ahead of 
the flame, a simple consideration of momentum balance at the flame leads 
to the result^ 


Pi = P 2 


(60) 


Hence, by the gas law, Pi/P2 ” equation ( 59 ) becomes 

“1 - - i) (6W 


But consideration of energy balance at the flame front shows that 


Q = Cp(T2 - Ti) 


(62) 


provided that — - « 1 . 

ai 


Hence, eqixation (6l) can be written simply as 


ui = 


5 t^ 

CpTi 


(65) 


Substituting equation (65) into equation (58) and making use of the gas 
law, it is found that the shock strength Pi/Pq must satisfy the 
equation 


7 - 1 

7 Pi/ Pi 


>a. f 

\Po J 


7+1 Pi ^ 7-1 
27 Po 27 


(64) 


5 por proof of eqs. (60) and (62), see, e.g., ref. 6 or 7 * 
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which after some simple reductions hecomes identical to equation (57) • 

Thus, it is seen that equation (57) indeed gives the correct value of 
the strength of the shock wave generated hy combustion. The fact that 
in this case the shock strength can be predicted exactly by another 
method has not made the usefulness of equation (59) less, for, in the 
first place, it has strengthened confidence in its application to 
practical problems, and, secondly, equation (59) can be applied to many 
other cases (some of which have been mentioned at the beginning of this 
section) where no other simple means is available for estimating the 
approximate strength of the shock wave produced. Finally, the corre- 
sponding problem in three dimensions can be solved by the theory 
developed in the last section. 

Mirels suggested that the condition under which the pressure wave 
produced by flame propagation can be found by considering an equal rate 
of heat release at a fixed station be considered. He noted that the 
equivalence of these two approaches is indicated for the planar case 
but is assumed for the three-dimensional case. (See the statement after 
equation (48): "If the flame speed is small, the heat may be thought 

of as being released at the center of the sphere.") 

In response to this suggestion, the author would like to add that 
before any comparison can be made for the three-dimensional case the 
exact solution for the flow field generated by a uniformly expanding 
spherical flame must be known. This solution can, in fact, be constructed. 
It may be worth while to point out that, depending on the values of the 
flame speed and the heating value of the mixture, the flow field may 
assume qualitatively different natures. When the flame speed and heating 
value of the mlxtiure are low enou^ (the present case), the flow field 
outside the flame front is similar to that generated by a uniformly 
expanding sphere, the speed of expansion of which is related to the flame 
speed and the heating value of the mixture. Inside the spherical flame 
the medium is at rest. When the flame speed and heating value of the 
mixture are hlg^ enou^, the flow field inside the spherical flame is no 
longer entirely at rest. There is now a family of central-expansion 
spherical waves following immediately behind the flame front which is 
now propagating at the lower Chapman- Jouguet speed. The flow field 
outside of the flame front is still similar to that generated by a 
uniformly expanding sphere. When the flame speed and heating value of 
the mixture are extremely high, the flame front catches up with the shock 
wave to form a detonation front ri^t after the mixture is ignited. This 
last case has been analyzed by Taylor (ref. 8) and lndei)endently by 
Doering and Bvirkhardt (ref. 9)- However, all these exact solutions 
cannot be given in closed analytic forms which involve only the elemen- 
tary transcendental functions. When the flame speed is small, the shock 
wave generated by the flame is extremely weak, mainly because of the 
fact that the shock wave is propagating into an open space. On the other 
hand, the equivalence of the flow field produced by the flame and a heat 
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source releasing heat at an equivalent parabolic rate is true only when 
the flame speed is small. Consequently, it does not seem worth while 
from the practical point of view to look into this particular example 
in any greater detail. 

When he was commenting on the equivalence of a flame front and a 
heater, the following argument was also advanced by Mirels: It can be 
shown for the planar case that the equivalence of the two approaches 
requires that the ratio of specific heats be the same for the burned 
and uribumed gases and that the kinetic-energy terms be negligible. Take 
the case of a flame originating at x = 0 at time t = 0. The tx-diagram 
is indicated in figure 11(a) . The equivalent problem, with heat addition 
at section x = 0, is Indicated in figure 11(b). Conditions in region (l) 
of the figures are the same for both cases but conditions in region (2) 
differ. In addition, the extent of region (2) is greater in figure 11(a) 
than in figure ll(b) since the flame moves faster than does the contact 
surface. However, the energy per unit volume of the gas behind the shock 
(for a given shock strength) can be shown to be the same for regions (l) 
and (2) and for figures ll(a) and 11(b) . Neglecting kinetic energy, the 
energy per unit volume is 


E = pC^T = 

7-1 

Since p Is constant behind the shock, the energy per unit volume is 
constant provided y is the same for the burned and imburned gases. 
Therefore the details of the temperature and density distributions 
behind the shock are unimportant. The strength of the shock depends on 
the rate of heat addition and the two approaches are equivalent. 

With regard to this discussion, the fact that a flame front and a 
heating element are dynamically equivalent^ when the ratio of specific 
heats is the same for the burned and unbiurned gases and when the flame 
speed is small compared with the local sound speed was Independently 
found by the author in his study of the mecheinism of generation of pres- 
sure waves at a flame front (ref, lO). Actually, the statement is rigor- 
ously true only if there is a current of flow through the heating element 
with velocity equal to the flame speed. The rigorous demonstration will 
not be presented here. 


^^en the pressure and velocity fields produced in two systems are 
identical, the systems are said to be "dynamically equivalent." Note 
that the temperature and density fields produced in the two systems need 
not be the same. 
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CONCLUDING REMARKS 


When the effect of heat conductivity is neglected, the flow field 
resulting from the addition of heat into a medixjin is caused by the 
volumetric expansion of the heated gas. When the rate of heat release 
is moderate, a "reduction theorem" can he derived which reduces the 
problem of heat addition in a single plane in a tube of constant cross 
section to a problem of piston motion in the tube (the plane of heat 
addition being perpendicular to the tube axis) . The fictitious pistons 
in the theorem correspond in reality to the interfaces which separate 
the heated and unheated gases. The solution for the more general case 
of heat addition in a region inside the tube can be constructed by 
superposition. 

The exact solution for the flow field produced by uniform heating 
in a plane at a constant rate is also given. In particular, the strength 
of the shock waves resulting from such heating is calculated in terms of 
the (constant) rate of heat release- The formula also gives the asymp- 
totic strength of the shock waves resulting from heating at a constant 
rate a finite volume of gaseous medium inside an infinitely long tube of, 
constant cross section. 

When the heat is added into the medium at a single point at a rate 
proportional to the time squared, the heated gas expands at a uniform 
rate, much like a uniformly expanding sphere. Taylor's solution enables 
one to calculate the relation between the shock wave produced and the 
rate of heat release- This relation also represents asymptotically the 
strength of a pressure wave generated by a closed flame front in a 
combustible mixture expanding uniformly at a constant speed which is 
small compared with the local sound speed. 


The Johns Hopkins University, 

Baltimore, Md., December 11, 


1955 ^ 
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Figure 1*-.- Shock strength versus parameter of rate of heat release 
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Q(t) 



Figure 5 .- Rate of heat release of a spark. 




(h) Velocity field. 

Figure 6 .- Pressure and velocity fields produced by a spark. 
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